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A method for obtaining a model of a surface, characterized by the steps which involve obtaining measurements of geomet- 
rical data concerning specific points on the surface, making a grid of the surface, with the grid passing through said points, mem- 
orizing, at an address which is specific to each node of the grid, the coordinates of the node, the number of satellites of the node, 
information for access to the addresses of said satellites and thereafter to information which relates to them, and geometrical data 
which may be associated with said node; for each node, defining a local roughness index obtained from a weighted sum of the 
current coordinates of the node and its satellites, defining the sum of an overall roughness index representing the sum of all the 
local roughness indices, and of an overall index of the infringement of said geometrical data, iteratively adjusting the coordinates 
of indefinite nodes, by using at each adjustment the sum of a weighted combination of current node neighbour coordinates and of 
a combination of geometrical data associated with the node, in order to minimize said sum, and creating a model of the surface 
on the basis of the adjusted coordinates. 

(57)Abrege 

L'invention concerne un procede d'obtention d'une moderation d'une surface. II est caracterise par les etapes consistant a 
obtenir des mesures de donnees geometriques relatives a des points specifiques de la surface, effectuer un maillage de la surface 
passant par lesdits points, memoriser a une adresse specifique a chaque noeud du maillage les coordonnees du noeud, le nombre 
de satellites du noeud, une information d'acces aux adresses desdits satellites et par suite aux informations y relatives et des don- 
nees geometriques eventuellement associees audit noeud; pour chaque noeud, definir un indice de rugosite local obtenu a partir 
d'une somme ponderee des coordonnees actuelles du noeud et de ses satellites, definir la somme d'un indice de rugosite global, 
constitue par l'addition de tous les indices de rugosite locaux, et d'un indice global de violation desdites donnees geometriques, 
ajuster de facon iterative les coordonnees des nceuds indetermines, en utilisant a chaque fois l'addition d'une combinaison ponde- 
ree des coordonnees actuelles de voisins du noeud et d'une combinaison des donnees geometriques associees au noeud, pour ame- 
ner ladite somme a un minimum, et creer une representation de la surface a partir des coordonnees ajustees. 
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1 

PROCEDE DE MODELISATION D'UNE SURFACE ET DISPOSITIF POUR 
SA MISE EN OEUVRE . 

La presente invention concerne d'une facon 
g6n6rale 1 r investigation dans des corps k trois 
dimensions, et concerne en particulier un nouveau proc6d6 
pour l'obtention d'une representation d'une surface 
situ6e au sein de ce corps k partir d'un ensemble limite 
de donn6es g6om6triques connues relatives k ladite 
surface, 

L 1 investigation dans des corps tridimensionnels 
est une preoccupation majeure tout particuliSrement dans 
les domaines de la g6ologie, de la gfeophysique et de la 
biologie. 

Ainsi , par exemple dans le domaine de la 
g6ophysique, on souhaite obtenir des representations 
aussi fiddles que possibles de surfaces situ£es par 
exemple k 1' interface entre deux zones de natures ou de 
propri6t6s diff6rentes, k partir de donn6es obtenues lors 
de campagnes de prospection, voire en cours 
d' exploitation d'une richesse souterraine. 

En effet, l'efficacite de 1 ' investigation du corps 
tridimensionnel , et en particulier de la recherche de 
p6trole ou autre richesse souterraine, est conditionn£e 
par la fid61ite avec laquelle on peut reconstituer et 
repr6senter ce genre de surface. 

Par ailleurs, dans le domaine medical et 
biologique, l'on sait actuellement obtenir, par des 
proc6d6s divers, des representations de coupes d'un corps 
vivant ou analogue, et l'on recherche egalement k 
reconstituer et k repr6senter, de fa<?on 
tridimensionnelle, une surface situ6e k 1' interface entre 
deux milieux, et notamment les contours d'un organe ou 
autre* 

II existe k l'heure actuelle un certain nombre de 
techniques de modeiisation permettant d' obtenir la 
representation d'une surface dans un corps k trois 
dimensions que l'on ch rche & exploiter ou k traiter. On 
citera en particulier les interpolations de B6zier ou les 
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fonctions splines (on se r6f6rera pour plus de details 
par exemple k 1'ouvrage Geometric' Modelling, M. E. 
Mortenson, Ed, John Wiley, 1985). 

Cependant, ces techniques connues sont toutes tr£s 
5 mal adapt6es k la prise en compte de donn^es h6t6rog£nes, 
k savoir des donn6es relatives k la fois aux coordonn6es 
de points de la surface, k des orientations impos6es de 
la surface, k des liens g6om6triques entre deux points 
distants, etc. En outre, ces techniques sont inadapt6es k 

10 

une trop grande dispersion des surfaces ou encore ct des 
anomalies de surfaces telles que des plis, des cassures, 
ou encore des surfaces discontinues. Plus pr6cis6ment, 
dans certains cas, les fonctions mises en oeuvre dans ces 
techniques soit ne convergent pas, soit n'admettent 

15 aucune solution, soit admettent au contraire une 
plurality de solutions. 

Enfin les techniques connues n % off rent pas la 
possibility de prendre en compte la notion de plus ou 
moins grande certitude avec laquelle certaines donn6es de 

20 la surface sont connues. 

La pr6sente invention vise k pallier les 
inconv6nients des proc6d6s de la technique ant6rieure et 
k proposer un proc6d6 qui puisse prendre en compte des 
informations g6om6triques non homog&nes et/ou trfes 

25 dispers6es, ainsi que les anomalies pr6cit6es et 
d'autres, en off rant k chaque fois une mod61isation 
unique et sans Equivoque de la surface. 

Un autre objet de 1* invention est de proposer un 
proc6d6 qui puisse s'appuyer sur des donnGes g6om6triques 

30 en prenant en compte simultan6ment un degr§ de certitude 
ou de precision li6 k ces donn6es. 

A cet effet, la pr6sente invention concerne tout 
d'abord un proc6d6 d'obtention d'une mod61isation d'une 
surface repr6sentative notamment de 1' interface entre 

35 deux zones de natures ou de propri6t6s diff6rentes dans 
un corps k trois dimensions tel qu'une. formation 
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g6ologique ou un corps vivant, caract6ris6 en ce qu*il 
comprend les 6tapes consistant k: 

- obtenir au moyen d' appareils de mesure un 
ensemble de donn6es g6om6triques relatives k la surface 
et associ6es respectivement k des points donn6s de ladite 
surface; 

~ effectuer un maillage de la surface de telle 
sorte que l f ensemble desdits points donn6s soit un sous- 
ensemble des noeuds du maillage; 

- m^moriser k une adresse donn6e d'une m6moire, 
sp6cifique k chaque noeud du maillage, les informations 
suivantes: 

* les coordonn6es du noeud consid6r6, 

* le nombre de noeuds satellites du noeud 
15 consider 6/ 

* une information permettant d*acc6der aux 
adresses sp&cifiques desdits noeuds satellites et par 
suite aux informations relatives k ces derniers, 

* le cas 6ch6ant des donn6es g6om6triques 
20 associ6es audit noeud consid6r6, 

- pour chaque noeud du maillage, d6finir un indice 
de rugosit6 local obtenu k partir d'une somme pond6r6e 
des coordonn6es actuelles du noeud et de ses satellites, 

- d6finir la somme d'un indice de rugosit6 global, 
25 constitu6 par 1" addition des indices de rugosit6 locaux 

associ6s k chaque noeud, et d'un indice global de 
violation desdites donn6es g6om6triques ; 

- ajuster de fagon it6rative les coordonn6es de 
chacun des noeuds dont on ne connalt pas les coordonn6es 

30 pr6cises, en utilisant k chaque fois I'addition d'une 
combinaison pond6r6e des coordonn6es actuelles des 
satellites et des satellites des satellites dudit noeud 
et d'une combinaison des donn6es g6om6triques associ6es 
audit noeud, de mani&re k amener ladite somme k un 

35 minimum; 



WO 91/04544 



4 



PCT/FR90/00676 



- cr6er vine repr6sentation de la surface k partir 
des coordonn6es ajust6es de chacun des noeuds. 

L* invention concerne 6galement un dispositif pour 
l'obtention d'une mod61isation d'une surface 
5 repr6sentative notamment de 1' interface entre deux zones 
de natures ou de propri6t6s diff6rentes dans un corps ci 
trois dimensions tel qu'une formation g6ologique ou un 
corps d'etre vivant, caract6ris6 en ce qu'il comprend : 

- des moyens de mesure un ensemble de donn6es 
10 g6om6triques relatives ci la surface et associ6es 

respectivement k des points donn6s de ladite surface; 

- des moyens pour effectuer un maillage de la 
surface de telle sorte que 1' ensemble desdits points 
donn6s soit un sous-ensemble des noeuds du maillage; 

15 - une m6moire contenant k des adresses donn6es, 

sp6cifiques k chaque noeud du maillage, les informations 
suivantes: 

* les coordonn6es du noeud consid6r6, 

* le nombre de noeuds satellites du noeud 
20 consid6r6, 

* une information permettant d'acc6der aux 
adresses sp6cifiques desdits noeuds satellites et par 
suite aux informations relatives k ces derniers, 

* le cas 6ch6ant des donn6es g6om6triques 
25 associ6es audit noeud consid6r6, 

- des moyens de calcul pour ajuster de fa<?on 
it6rative les coordonn6es de chacun des noeuds dont on ne 
connait pas les coordonnfees pr6cises, en utilisant k 
chaque fois 1' addition d'une combinaison pond6r6e des 

30 coordonn6es actuelles des satellites et des satellites 
des satellites dudit noeud et d'une combinaison des 
donnfees g6om6triques associfees audit noeud, de manidre k 
amener k un minimum une somme d'un indice de rugosit<§ 
global, constitu6 par 1' addition d* indices de rugosit6 

^ locaux respectivement associ6s k chacun des noeuds et 
chacun obt nu k partir d'un somme pond6r6e des 
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coordonn6es actuelles dudit noeud et de ses satellites, 
et d'un indice global de violation desdites donn£es 
g&om&triques ; et 

- des moyens pour cr£er une representation de la 
5 surface h partir des coordonn6es ajust6es de chacun des 
noeuds . 

D'autres aspects , buts et avantages de la prfesente 
invention apparaitront mieux & la lecture de la 
description d£taill£e suivante d'une forme de realisation 
10 pr6f6r6e de celle-ci, donn6e k titre d' exemple non 
limitatif et faite en reference aux dessins annexes, sur 
lesquels : 

la figure 1 est un exemple de graphe incorporant 
un maillage utilise pour mod6liser une surface, ce graphe 

15 etant d£fini par 1" ensemble de tous les c6t£s de facettes 
dont les sommets constituent 1* ensemble Q d£crit plus 
loin, et illustre 6galement des filches symbolisant des 
d6placements de noeuds du maillage effectu£s 
interactivement par 1 'utilisateur pour fagonner la 

2o surface; 

les figure 2(a) k 2(d) illustrent quatre cas de 
contr61e de la forme d'une surface d l'aide de quatre 
types diff6rents de contraintes vectorielles; les 
vecteurs sont supposes connus directement sur les 

25 figures 2(a) & 2(c) et connus par leur orthogonality & un 
vecteur V donn6 sur la figure 2(d); 

la figure 3 illustre un exemple de surface 
g£ologique complexe (d6me de sel) modeiis6e & l'aide du 
proc6d6 de la prfesente invention; 

la figure 4 illustre un exemple de surface 
biologique complexe (cerveau d'embryon) mod£lis6e & 
l'aide du proc6d£ de la prfesente invention h partir de 
donn6es ponctuelles obtenues par une succession de 
coupes ; 

la figure 5 illustre 1 'organisation de la m£moire 
associ6e h un atome A(k) de noyau <f>k ; 
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la figure 6 illustre un premier exemple de 
memorisation d'une orbite d'un atome donne, dans un 
tableau de nn cases? 

la figure 7 illustre un second exemple de 
5 memorisation d'une orbite d'un atome donne, dans nu 
tableaux chain6s de deux cases chacun; 

la figure 8 illustre la memorisation d'une donn6e 
gfeometrique ou contrainte; 

la figure 9 illustre la mani6re de m6moriser un 

10 

ensemble de donn6es g6om6triques ou contraintes relatives 
k un atome A(k) donn6? 

la figure 10 illustre la memorisation dans un 
tableau de 1' ensemble des atomes relatifs A une surface 
S; et 

15 la figure 11 illustre 1 ' enchalnement de l'acces 

audit ensemble d' atomes. 

On indiquera tout d'abord que, d'une figure & 
l 1 autre, des 616ments ou parties identiques ou similaires 
sont d6sign6s par les mSmes num6ros de reference. 

20 En outre, on va indiquer ci-dessous un certain 

nombre de publications dont les contenus respectifs sont 
4 consid6rer comme inclus dans la pr6sente description. 

(1) A method of Bivariate interpolation and smooth 
surface fitting based on local procedures, H. AKIMA, J. 

25 ACM 17.1, 1974? 

(2) Shape reconstruction from planar cross-sec- 
tions, J.D. BOISSONNAT, ACM Transactions on Graphics, 
vol. No. 2, 1986; 

(3) Machine contouring using minimum curvature, 
30 I.C. BRIGGS, Geophysics 17,1, 1974; 

(4) Triangular Bernstein-B6zier patches, G. FARIN, 
Computer-aided geometric design, 1986; 

(5) Primitives for the manipulation of general 
subdivisions and the computation of Voronoi diagrams, L. 

35 GUIBAS et J. STOLFI, ACM Transactions on Graphis, vol. 
No. 2, 1985; 
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(6) Three-dimensional graphic display of disconnected 
bodies, J.L. MALLET, Mathematical Geology, vol. 20, no. 8, 
1988; 

(7) Geometrical modelling and geostatistics , J.L. 
5 MALLET, Proceedings of the Third International 
Geostatistics Congress/ Kluwer Academic Publishers, 1989? 

(8) Discrete Smooth Interpolation, J.L. MALLET, ACM 
Transactions on Graphics, April issue, 1989? 

(9) Geometric Modelling, M.E. MORTENSON, John Wiley, 
10 1989. 

Par ailleurs, la description qui suit d'une forme de 
realisation pr6f6r6e de 1' invention est elle-mdme suivie 
d'une annexe exposant une m6thode de reconstitution d'une 
surface d6riv6e de la methode expos6e dans le document (8) 
15 et sur laquelle est bas6 le proc6d6 de la pr6sente 
invention. 

1. DESCRIPTION DU PROCEDE DE CODAGE ET D 1 AJUSTEMENT 

DE SURFACES COMPLEXES SELON L* INVENTION 
1.1. Introduction 
20 1.1.1 Ob jet de la description 

On va d6crire ci-dessous un proc6d6 industriel de 
codage et d'ajustement de surfaces complexes rencontr6es 
par exemple lors de la modfelisation d' interfaces entre 
zones de natures ou de propri6t6s diff6rentes de couches 
25 g6ologiques ou de corps d'organismes biologiques. 

Le r6sultat de ce codage peut §tre transcrit sous forme 
d' images et/ou de maquettes en materiel solide (par exemple 
plastique) & l'aide d'un ordinateur graphique ou de tout 
autre proc6d6. En g6ologie, par exemple, ces images et/ou 
30 maquettes sont ensuite utilis6es pour optimiser 
1 'exploration , 1 'exploitation et la gestion des ressources 
du sous-sol telles que : 
. p6trole 
. minerai 
35 . eau 

. etc. 
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Le proc6d6 d6crit peut Stre decompose en deux sous- 
proced^s complementaires : 

1. un sous-procede de codage d'une surface 

2. un sous-procede d'ajustement d'une surface aux donnees disponibles. 
5 En geologie, par exemple, ces donnees peuvent etre : 

• la position (exacte ou approchee) de points sur une coupe de la 
surface pax un plan 

• la position (exacte ou approchee) de Intersection d'un forage avec 
la surface geoiogique correspondant a ia separation entre deux couches 

10 • le plan tangent (exact ou approche) en un point donne de la surface 

correspondant a la separation entre deux couches 

• le rejet (exact ou approche) d'une faille (=cassure) affectant la 
surface correspondant a la separation entre deux couches 

• etc ... 

1 5 Le sous-procede de codage est con$u de fa$on a optiiniser ies performances 
du procede d'ajustement base sur une variante de ia methode DSI (cf ref. 
[8]) specialement adaptee a la modeiisation geometrique des surfaces. Cette 
variante de la methode DSI utilisee ici n'a encore jamais ete publiee et est 
presentee en annexe. 

20 1.1.2 Decomposition d'une surface S en facettes 

Soit S une surface a modeliser. Comme le montre la Figure 4 et ainsi qu'il 
a ete suggere dans les references [6] et [7], celle-ci sera supposee approchee 
par un ensemble de facettes triangulaires et/ou polygonales. Ces facettes sont 
caracterisees par : 

25 • La connaissance des coordonnees de leurs sommets constituant un en- 
semble fini de N points . . ,$ N } regulierement repartis sur 5. Ces 
sommets sont numerates de 1 jusqu'a N et ies coordonnees du fc eme 
sommet $ k seront notees (y?f 

• La donnee de leurs cotes constitues par des liens entre couples de som- 
30 mets 

L'ensemble des cotes des triangles et/ou poiygones composant une surface 
S constitue un graphe G dont les noeuds sont precisement les sommets des 
triangles et/ou poiygones. Pour cette raison, et afin d'eviter des periphrases 

FEUILLE DE ^EMPLACEMENT 
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inu tiles, nous utiliserons de fa$on equivaiente le mot "noeud" a la place de 
"sommet de triangle et/ou polygone" : 

noeud = { sommet de triangle et/ou polygone } 
1.1.3 Definitions preliminaires 

5 Ensemble des satellites A(fc) 

Soit (p k le fc eme noeud constituant une surface 5. Nous appellerons ensemble 
des satellites de <p k et nous noterons A(fc) 1'ensemble des noeuds tpj differents 
de (p k tels que (£y f $ k ) soit un cote de triangle ou de polygone. En d'autres 
termes, 1'ensemble A(fc) est tei que : 

10 1) A(fc) = {^ il> ...,(? infc } 

(Vj e ^ $k sont * es extremites d'un 
meme cote de triangle ou de poly- 
gone 

Comme indique dans la relation ci-dessus, nous conviendrons de noter n k le 
15 nombre de satellites rattaches au noeud <p k . 

Cette notion d'orbite introduite dans les articles [6] et [7] n'a jusqu'a 
present pas fait l'objet de procede de codage particulier; une partie impor- 
tante de 1' invention consiste k proposer un proc6d6 de 

codage pour /i(k). 

20 Voisinage N(k) 

Soit <p k le fc eme noeud d'une surface S. Nous appellerons voisinage de (p k et 
nous noterons N(k) 1'ensemble des noeuds constitues par la reunion de $ k et 
de ses satellites. En d'autres termes, on a : 

N(k) = {&} U A(fc) 

25 Nous appellerons voisin du noeud $ k tout noeud appartenant a Tensembie 
N(k). En d'autres termes, on a : 

€ N(k) <=> = Voisin de <p k } 

Cette notion de voisinage introduite dans les articles [7] et [8] n'a jusqu'a 
present pas fait l'objet de procede de codage particulier; une partie importante 
30 de 1' invention consiste k proposer un procede de cndage 
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Atome A(k) et noyau associe 

Soit (p k le k eme noeud d'une surface S. Nous appellerons atome de noyau 
<p k et nous noterons A(k) le couple A*(fc)) constitue par la juxtaposition 
(et non la reunion) : 

• du noeud <p k 

• de l'ensemble A*(fc) des adresses ou des codes {A*(ji), . . . , A*(j„J} per- 
mettant de retrouver dans une memoire Pensemble des atomes {A(ji), 
A (j nk )}correspondants dont les noyaux sont les satellites de tp k 

En d'autres termes, on a : 

m = (^.A*(t))-(ft,{A*ci.),...,A*cjo» 

rik = nombre de satellites de <p k 

avec < A*(j) = Adresse ou code d'acces de i'atome 

A(j) correspondant au noyau du j eme 
satellite de (p k 

1 5 Cette notion d'atome introduite dans Particle [6] n'a jusqu'a present pas fait 
l'objet de procede de codage particulier; une partie import ante de 
1' invention consiste h proposer un proc§d6 de codage pour k(k). 

Contraintes et types de contraintes 

Nous appellerons contrainte relative a un atome A(fc), toute condition relative 
20 a la position du noeud <j> k correspondant au noyau de cet atome. Par exempie, 
en geologie on peut citer les contraintes suivantes qui sont d'une importance 
capitale pour Tajustement d'une surface aux donnees sures ou approchees 
observees ou mesurees par des precedes manuels ou geophysiques : 

• Les coordonnees (<p k ><p y ky du noeud <p k sont connues de fa$on exacte. 
25 Une telle contrainte est dite de type noeud de controle et est utilisee 

en geologie par exempie pour coder ^intersection de la surface modelisee 
avec un forage. 

• Les coordonnees fcfi du noeud tp k sont connues de fa$on ap- 
prochee avec un coefficient de certitude donne. Une telle contrainte est 

30 dite de type noeud de controle flou et est utilisee pour coder la posi- 

tion approximative d'un noeud de la surface modelisee supposee mesuree 
de fagon approximative par un procede donne. Par exempie, 
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- en geologie, les donnees sismiques sont de ce type, 

— en biologie, les donnees de coupes d'organes observees sous micro- 
scope ou par echographie ou encore par scanner sont de ce type. 

# Le vecteur {(p x — $k) joignant le noeud tp k a un autre noeud (p x est 
5 connu de fa$on approchee avec un coefficient de certitude donne. Une 

telle contrainte est dite de type lien vectoriel flou et est utilisee en 
geologie par exemple pour coder le rejet d'une faille, le noeud ^etant 
sur un bord de la faille et le noeud $ x etant aiors sur le bord oppose 
(voir Figure 2). 

10 • Le vecteur V orthogonal au plan tangent a la surface S au point corre- 
spondant au noeud <p k est connu de fagon approchee avec un coefficient 
de certitude donne. Une telle contrainte est dite de type normale floue 
et est utilisee en geologie par exemple pour coder le pendage des couches 
mesure par des procedes manuels ou geophysiques (pendagemetrie). 

15 # etc ... 

Les coefficients de certitude mentionnes ci-dessus sont des nombres positifs 
dont la valeur est censee representer la confiance que Ton peut avoir "a priori" 
dans les donnees qui lui sont attachees. La prise en compte des contraintes 
floues pr6sent6es ci-dessus constitue I'une des 
20 caractferistiques de 1' invention est pr6sent6e en detail 
dans 1' annexe (voir annexe - section 2.5 et 2.7.6.). 

Convention 

Afin de simplifier les notations, on convient d'identifier les numeros des noeuds 
d'une surface S avec les noeuds eux-memes. II s'ensuit que Ton a : 

25 « A(fc) = ensemble des satellites {$j x ,- . - i$j nk } ^ e Vk ou ensemble des 
numeros {j u . . . , j n|t } correspondants 

* N(k) = ensemble des voisins {^ 9 ^ & » • • • ? ^ e $k ou ensemble des 
numeros {A;, ji, . . . , j„ k } correspondants 

1.1.4 Etat de Tart 

30 Les methodes classiques 

On trouvera dans la reference [7] une description tres complete de Tetat de Tart 
concernant la modelisation de surfaces complexes. En particulier, dans cet 
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article, on expose les raisons pour lesqueiles les methodes classiques utilisees 
en cartographie automatique et en Conception Assistee par Ordinateur (C AO) 
ne sont pas adaptees a la modeiisation de telles surfaces comme on aurait 
pu le croire a priori. A la fin de cet article [7], on introduit brievement les 
5 fondements de la methode DSI (quelquefois appeiee DSI/DSA) permettant de 
surmonter les difficuites inherantes aux methodes classiques. 

La methode DSI 

Cette methode est basee sur la notion de "rugosite" d'une surface S composee 
de facettes triangulaires et/ ou poiygonales. Tout d'abord, la notion de rugosite 
1 0 R{<p\k) au noeud (p k d'une surface S est definie a Taide de la formule suivante 
(Voir ref. [7], [8] et Annexe) : 

R(m hi E a*)- a. ii a 

a€AT(fc) 

Les coefficient positifs ou nuis {u 0 (fc)} intervenant dans la formule ci-dessus 
sont des ponderateurs choisis par Tutilisateur charge de modeliser la sur- 
1 5 face. Parmi Tinfinite de choix possibles, Tun des plus interessants appeile 
"ponderation harmonique" consiste a poser : 



= | -|A(*)| si a = k 



si a € A(fc) 

Les rugosites locales ainsi definies sont ensuite combinees en une rugosite 
20 globale R(<f) definie comme suit : 

W)«Ert*)-WI*) 

fc 

Les coefficients positifs ou nuls /i(fc) intervenant dans la formule ci-dessus 
servent a moduler la contribution des rugosites locales R($\k) a la rugosite 
globale i?(<£). Si i'on decide par exempie d'utiliser une ponderation uniforme, 
25 il suffit de poser : 

/z(Jb) = l Vfc 

Parmi les autres choix possibles pour les coefficients fi{k) y on peut decider 
d'utiliser par exempie les ponderateurs suivants ou m est une constante donnee : 

/r x _ j 1 si (p k est un point de donnee 

30 "~ \ ra > 1 si <p k n'est pas un point de donnee 
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Le principe de la methode DSI consiste 2i ajuster la position 
de chacun des noeuds 9* de la surface S Ood61is6e de telle 
sorte que celle-ci soit la plus lisse possible tout en res- 
pectant les donndes et contraintes rggissant sa forme. 
5 L 1 article [8] mentionnd en bibliographie est consacrd k une 
presentation mathdraatique approfondie de la methode DSI et 
propose une technique de minimisation de la rugositd globale 
respectant les <donn£es et contraintes. 

Version interactive de la methode DSI 

1 0 La methode de minimisation de la rugosite globale de la methode DSI proposee 
dans Particle [8] est mal adaptee a une utilisation interactive sur calculateur 
electronique (ordinateur ou microprocesseur specialise), c'est pourquoi nous 
proposons en annexe une autre technique de minimisation originale sur laque- 
lie est basde 1 1 invention . Ainsi qu'il a ete dit precedemment, Tobjet de 

1 5 Invention est precisement de proposer un procede de codage et un procede 
d'ajustement de la surface S suivant cette nouveile version de la methode DSI. 

1.2 Exemples damages d'objets 

codes et ajustes par notre procede 

1*2.1 Exemple d 'image de surface geologique 
20 obtenue par notre procede 

La recherche du petrole est generalement effectuee de fagon indirecte en ce sens 
que plutot que de rechercher le petrole, on cherche des couches geologiques 
reservoir suceptibles de le contenir ou de le capturer; par exemple, la deformation 
des couches de sei donne naissance a ce que Ton appeile des "domes de sel" qui 

25 sont consideres comme generant d'excellents pieges a petrole. Pour pouvoir 
localiser avec precision la position du petrole il est necessaire de connaitre 
la forme exacte des couches "pieges" et malheureusement dans de nombreux 
cas il n'existait jusqu'a ce jour aucun procede pratique permettant de les 
modeliser de fagon satisfaisante; le cas des domes de sel est a ce titre considere 

30 comme Tun des plus complexes et pour prouver l'efficacite de notre procede 
de codage et ajustement nous presentons sur la Figure 3 Timage d'une telle 
surface obtenue avec celui-ci. 

FEUILLE DE REMPLACEMENT 



WO 91/04544 



- 14 - 



PCT/FR90/00676 



1.2.2 Exemple d 'image de surface biologique 
obtenue par notre precede 

Les procedes d'acquisition de donnees tels que scanner ou echographie utilises 
en medecine sont mathematiquement identiques aux procedes d'acquisition de 

5 donnees sismiques utilises en prospection petroliere. Par exemple, Techographie 
est identique a la sismique reflexion et le fonctionnement des scanners est 
base sur la tomographic tout comme la tomographic sismique. Des lors, il 
n'est pas etonnant que notre procede de codage et ajustement des surfaces 
soit egalement utiiisable pour modeiiser la peau d'organes biologiques et nous 

1 0 proposons sur la Figure 4 le modele du cerveau d'un embryon humain de 7 mil- 
limetres de long obtenu a Taide de notre procede de codage et d'ajustement. 
II est k noter qu'un aussi petit organe (le cerveau fait 
1 millimetre de long) ne pouvait jusqu'<k present §tre 
observe que sous forme de coupes microscopiques . 



1.3 Procede de codage d'une surface 5 

1.3.1 Notion de memoire et d'adresse de memoir e 

Nous appellerons "memoire* tout dispositif de calculatrice electronique qui 
enregistre des donnees, des codes ou resultats partieis qu'il restitue au mo- 
20 ment opportun pour les faire intervenir dans la suite des operations. De plus, 
nous appellerons "adresse" d'une memoire, tout code permettant de retrouver 
l'emplacement physique de la memoire a Tinterieur de Tordinateur. 

1.3.2 Codage d'un atome A(k) 

Nous proposons de coder et stocker chaque atome A(k) dans un ensemble de 
2 5 memoires a l'aide du procede suivant ou Tune de ses variantes : 

Procede de codage 

Comme le suggere la Figure 5 , dans une zone memoire contigue commenqant 
a Tadresse on stocke de fa?on consecutive les informations suivantes 

relatives a Patome A(k) correspondant au noeud $ k : 

30 1. n yau = 
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# soit les coordonnees (v?f iViiVk) noeud <p k 

# soit l'adresse ou le code ou toute autre information permettant de 
retrouver les coordonnees (v^V'/kv 5 *) ^ u noeud (p k 

2. nbjsat = 

5 # soit le noinbre n* de satellites lies du noeud tp k 

# soit l'adresse ou le code ou toute autre information permettant de 
retrouver le nombre n k de satellites lies du noeud (p k 

3. sat = 

Adresse d'une zone memoire contenant les adresses ou les codes d'acces 
1 0 ),..., A*{j nu )} permettant de retrouver les n* atomes {A(ji ),..., 

A 0nJ} dont ies n °y aux {?j x i • • • yfjnj sont les satellites du noeud <p k 

4. const = 

Codage de 1 'ensemble des contraintes (cf section 1.1.3) pouvant afFecter 
le noeud Ce codage sera presente en details dans la section 1.3.3. 

15 5. info = 

Zone memoire de taille variable contenant zero, une ou plusieurs infor- 
mations complementaires concernant l'atome A(k). Ces informations 
sont specifiques a 1'application particuliere envisagee; par exemple en 
prospection petroliere, si A(k) est un atome correspondant a la surface 
20 separant deux couches C\ et Ci, alors il peut etre interessant de stocker 
dans le champ info des informations telles que : 

# liste des proprietes physiques (porosites, vitesse des ondes sismiques, 
...) de la couche C\ au noeud <p k correspondant a l'atome A(k) 

# liste des proprietes physiques (porosites, vitesse des ondes sismiques, 
25 ...) de la couche au noeud $ k correspondant a l'atome A{k) 

# liste des proprietes geologiques (fades geologique, presence de petroie, 
...) de la couche C\ au noeud (p k correspondant a l'atome A(k) 

# liste des proprietes geologiques (fades geologique, presence de petroie, 
...) de la couche C 2 au noeud (p k correspondant a l'atome A(k) 

30 * etc 
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Variante 1 

La zone memoire dans laquelle sont stockees les adresses ou les codes d'acces 
{A*(ji), . . . , A*(j nk )} permettant de retrouver les n k atomes {A(j x ), . . . , A(j n J} 
dont les noyaux sont les satellites {c^, . . . , £j nfc } du noeud <p k peut etre struc- 
5 turee de deux fagons : 

1. La premiere fa$on illustree par la Figure 6 consiste a utiliser un tableau 
de nfc memoires consecutives contenant les adresses {A*(ji), . . . , A^jnJ} 

2. La seconde fa$on illustree par la Figure 7 consiste a utiliser tableaux 
comportant chacun 2 memoires consecutives telles que pour le a eme de 

1 o ces tableaux : 

# la premiere memoire contienne Tadresse A*(j a ) de Tatome corre- 
spondant au a eme satellite de l'atome A(k) 

# la seconde memoire contienne I'adresse du (a+l) eme tableau. Si a 
est egal au nombre de satellites (nb«sat=njt), alors cette seconde 

1 5 memoire est soit inutilisee, soit contient un code indiquant que Ton 

a atteint ie dernier satellite de A(k) 

Variante 2 

II est possible de ne pats stocker le champ nb jsat mais dans ce cas, si la liste des 
adresses ou codes d'acces {A*(j\) y . . . , A*(j nfc )} est codee sous forme d'un seul 
20 tableau (voir variante 1), il faut rajouter un emplacement memoire A*(n k + 1) 
dans lequel on place un code indiquant que la liste A*(fc) est terminee. 

Variante 3 

Le champ sat peut etre choisi suffisamment grand de fa§on a y stocker di- 
rectement les adresses ou codes d'acces {A*(jt), . . . , A*(j nk )} 

25 Variante 4 

Un grand nombre de variantes equivalentes sont possibles suivant la taiile 
du champ info et de la fagon dont on redecoupe celui-ci. II est toutefois 
bien entendu que ces variantes n'ont aucune importance pour l'utiiisation 
de l'atome A(fc) en connection avec ralgorithme .d'ajustement base sur la 
30methode DSI decrite en annexe. 
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Variante 5 

On peut rajouter d'autres champs a I'atome A(k) decrit ci-dessus mais ceux-ci 
n'ont aucune importance pour l'utilisation de i'atome A(k) en connection avec 
Talgorithme d'ajustement base sur la methode DSI decrit en annexe. 

5 Variante 6 

Quelles que soient les variantes, la seule chose importante pour Hmplementation 
efficace de Palgorithme DSI est la suivante : 

Ayant I'adresse ou Ie code d'acces A*(k) d'un atome A(k), il est 
necessaire d'acceder de la fa^on la plus directe possible aux atomes 
10 dont les noyaux sont les satellites du noyau de A(k). 

I/object essentiel du procede de codage de A(k) dans une memoire decrit 
ci-dessus ainsi que toutes ses variantes est precisement de permettre Tacces 
direct aux satellites. 

1.3.3 Codage des contraintes 

1 5 Chaque contrainte relative a un atome A(k) sera stockee dans une zone memoire 
specifique dont la description est proposee dans ce qui suit : 

Procede de stockage d'une contrainte 

Le codage d'une contrainte particuliere relative a un atome donne est effectue 
dans une zone memoire comprenant les deux champs contrainte et next 
20 comme Pindique la Figure 8 . Ces deux champs sont utilises de la fa?on suiv- 
ante : 

« contrainte = 

Zone memoire contenant : 

- soit les informations relatives a la contrainte codee (voir section 
25 plus bas) 

— soit Tadresse ou le code d'acces d'une zone memoire contenant ef- 
fectivement les informations relatives a la contrainte codee 

Parmi les informations relatives a la contrainte codee, doit toujours fig- 
urer un code permettant de reconnaitre le type de la contrainte codee 
30 (cf section 1.1.3). 
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• next = 

zone memoire contenant : 

- soit un code indiquant que la contrainte consideree est la derniere 
des contraintes affectant l'atome A(k). 

5 - soit l'adresse d'une zone memoire contenant ia contrainte suivante 

(voir section suivante). 

Procede de stockage de toutes les contraintes affectant un atome 

Le codage d'un atome A(k) decrit dans la section 1.3.2 contient un champ 
const que nous utiliserons pour stocker : 

1 0 • soit un code indiquant que l'atome A(k) n'est affecte par aucune con- 
trainte 

• soit l'adresse d'une zone memoire contenant la premiere contrainte af- 
fectant l'atome A(k) 

Si l'adresse de la n eme contrainte relative a l'atome A(k) est stockee dans le 
15 champ next de la (n — l) eme contrainte (n > 1), alors, comme le suggere la 
Figure 9, on est assure que le procede de codage decrit permet detacher a 
l'atome A(k) toutes les contraintes qui lui sont relatives. Ce mode de codage 
permet de raj outer ou de supprimer faciiement des contraintes sans remettre 
en question l'organisation des autres donn&s. 

20 Exemples d'informations relatives a une contrainte donnee 

Dans la section 1.1.3 nous avons presente a titre d'exemple trois types de con- 
traintes tres importants pour le codage et Pajusternent des surfaces complexes 
rencontrees en biologie ou geoiogie. A chacun de ces types de contraintes cor- 
respondent des informations qui doivent etre associees au champ contrainte 
25 decrit au paragraphe precedent et schematise sur la Figure 9: 

• Pour une contrainte de type noeud de controle, en plus du type, il est 
necessaire de stocker au minimum les donnees suivantes : 

— les coordonnees (ip% y v^^fc) approchees du noeud (p k correspondant 
au noyau de l'atome A(k) auquel est attachee la contrainte 

30 — le coefficient de certitude (nombre positif) mesurant ia confiance 

que Ton peut avoir dans les coordonnees (^^V 5 !)- 
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# Pour une contrainte de type lien vectoriel, en plus du type, il est 
necessaire de stocker au minimum les donnees suivantes : 

- i'adresse du noeud <£ A auquel est lie le noeud (p k correspondant au 
noyau de l'atome A(k) auquel est attachee la contrainte 

5 — les trois composantes approchees du vecteur ((p x — joignant le 

noeud (p k a au noeud tp^ 

- le coefficient de certitude (nombre positif ) mesurant la confiance 
que Von peut avoir dans les composantes approximatives du vecteur 

(ft - ft) • 

1 0 # Pour une contrainte de type normale, en plus du type, il est necessaire 
de stocker au minimum les donnees suivantes : 

- les trois composantes approchees du vecteur V orthogonal a la sur- 
face 5 au noeud $ k correspondant au noyau de l'atome A(k) auquel 
est attachee la contrainte 

15 - le coefficient de certitude (nombre positif) mesurant la confiance 

que Von peut avoir dans les composantes approximatives du vecteur 
V. 

La presentation detaillee de ces trois types de contraintes est proposee dans 
les sections 2.5 et 2.7.6 de ['annexe. 

20 Variant e 

Bien que Ton n'ait parle dans ce qui precede que d'un seui coefficient de 
certitude pour les composantes en x } y et z de chaque contrainte, on pourra 
affecter a ces 3 composantes 3 coefficients de certitude differents si necessaire. 

1*3.4 Codage d'une surface 5 

2 5 Afin d'optimiser la programmation de la variante de la methode DSI presentee 
en annexe, nous proposons de coder une surface 5 composee de facettes tri- 
angulares et/ou polygonales sous forme d'un ensemble de N atomes dont 
les noyaux (noeuds) sont les sommets de ces triangles et/ou polygones. Les 
donnees et codes correspondant a chacun de ces atomes {A(l ),..., A( AT)} 

30 peuvent etre stockes : 

• soit dans un tableau corame le suggere la Figure 10 
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# soit dans le champ at m de 2ones memoires chainees entre elles par leur 
champ next comme le suggere la Figure nJDans ce cas, le champ next 
associe a Tatome numero k doit contenir Tadresse ou le code d'acces * 
de Tatome suivant numerote (k + 1); le champ next associe au dernier 

5 atome A(N) doit contenir un code indiquant que le dernier atome a ete 

atteint. 

La deuxieme solution est la plus interessante si Ton souhaite pouvoir aj outer 
ou supprimer facilement des atomes, par contre elie consomme d'avantage de 
place memoire que la premiere. 

10 1*3.5 Variantes 

Pour des commodites de presentation, les zones memoires utilisees pour coder 
des informations ont ete decoupees en champs et chacun de ces champs a ete 
attribue : 

• un nom propre 

15 # une place dans la zone memoire correspondante 

II est bien evident que le procede de codage decrit dans ce brevet ne depend 
pas des noms ni de Tordre de ces champs a Tinterieur des zones memoire 
utilisees pour le codage. Par ailleurs, il est toujours possible de rajouter des 
champs a des fins specifiques a certaines applications sans que cela modifie le 
20 procede de stockage. 

1.4 Ajustement d'une surface 5 

1.4.1 Introduction 

Etant donne une surface 5 codee suivant le procede decrit dans la section 
1.3, nous nous proposons dans ce qui suit de decrire un procede base sur ce 
2 5 codage et permettant d'ajuster de fagon optimale les coordonnees de certains 
des noeuds . . • >^iv} et cfici ^ e sorte Q ue ^ a surface 5 : 

• soit la plus lisse possible au sens du critere de rugosite globale de la 
methode DSI (voir references [7], [8] et annexe ainsi que la section 1.1.4) 

# respecte le mieux possible les informations et contraintes connues avec 
3 0 plus ou moins de certitude concernant la forme de cette surface ( voir 

sections 1.1.1 et 1,1.3) 
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15 



20 



Pour ce faire, on utiiisera une variante originale de la methode DSI presentee 
en annexe specialement congue pour utiiiser ie codage decrit dans la sec- 
tion 1.3. Compte tenu des informations disponibles observees ou mesurees, 
cette methode ajuste successivement et de fagon iterative les trois coordonnees 
(<r>t,<Pa><pl) de chaque noeud <£ a non fixe suivant une formule du type : 



nouvelle valeur de ip% 
nouvelle valeur de \p\ 
nouvelle valeur de y>* 



/*( anciennes valeur s de <p* , . . . , <p%) 
fl{ anciennes valeurs de <p\, , . . , tp y N ) 
/*( anciennes valeurs de <p\^ . . . , tp* N ) 



Plus precisement, si Ton se reporte a I'annexe presentant la variante utilisee de 
1 0 la methode DSI, ces fonctions d'ajustement /£(), f%() et /*() sont deduites de 
la forme locale des equations DSI au point numero a et ont la forme suivante 
(voir sections 2.4.1, 2.4.2 et 2.7.3 de I'annexe) : 



£(•••) = 



l 



MS 



{ 



E {/<(*) ""W • E »'(*)• 9l>} 



E{rr - Qr)} 



i 

£(•••) = -k- { E {*(*)«•(*)• E AQ'rt) 

+ E{rr-gn} 



avec 



25 



f Mr = M a +E7f a 

Mr=M a +E7r 
M a = e M*) (*•(*))' 

fcew(a) 
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15 



La nature et le role des coefficients de ponderation {u a (fc)} et fi(k) ont ete 
indiques section 1.1.4 et sont detailles dans I'article [8] ainsi qu'en annexe. Les 
valeurs des termes {Tf, If, Tf a }, {yf a , y? a , 7?"} et des termes {Qf a , QF, Q?} 
dependent des types et des valeurs des contraintes affectant le noeud (p a et 
leur formulation exacte est presentee en annexe. 

Dans le cas particulier ou les coefficients de ponderation {v a (k)} corre- 
spondent a la ponderation harmonique presentee section 1.1.4 et en annexe, 
les fonctions de ponderation /*(), /£() et /*() se simplifient et prennent la 
forme suivante : 



20 



25 



0 =J_.f 



M«) ' |A(a)| • ( E Vf) 

V fc6A(o) 7 

E M*)-{(E v5)-|A(t)|-vs] 

fc6A(er) l V /S6A(*) 7 } 

EUr-GT) ] 



£(•••) = 



M*' 



[ 



ti«) • |A(a)| -(Erf) 

- E **(*)•{( E vS) - 1A(fc)| - vs] 

*€A(a) 1 \»€A(») 7 J 

- sw-on ! 

fi(«)'|A(a)|.( E rf) 

X fc6A(o) 7 

- E M*0-{(E V*)-|A(*)|-rf} 

fc€A(a) 1 X <J6A(*> 7 "* 

- Effr-tf"} ] 



avec : < 



f Mr = M a +E7f 
M^ = M a + E7? 

M. - Ma) • |AW + E M*) 

*€A(o) 



2a 
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Nous proposons dans ce qui suit un procede de calcul de ces fonctions base 
sur le codage presente au paragraphe 1.3: 

1.4.2 Notations 

Afin de simplifier la description des procedes de calculs presentes dans ce qui 
5 suit, il est commode d'introduire les notations et definitions suivantes : 

• Par definition, nous appellerons bloc d'operations, toute liste cooperations 
precedee d'une accolade ouvrante w {" et terminee par une accolade fer- 
mante "}". Les operations figurant dans un bloc d'operations peuvent 
etre des operations element aires et/ou d'autres blocs d'operations. 

1 o • Par definition, nous appellerons memoire de travail toute memoire utilisee 
pour stocker des resultats intermediates au cours d'un procede de calcul. 

• Par definition, nous appellerons '"chargement" de Texpression e dans 
une memoire m ^operation consistant a evaluer la vaieur de e puis a 
ranger le resuitat dans la memoire m. Par convention, une telle operation 

1 5 sera notee de la fagon suivante : 

m «~ e 

• Pour toute memoire m utilisee pour stocker une vaieur numerique, nous 
appellerons "incrementation" de m par 1'increment i Toperation consis- 
tant a calculer la vaieur m + i puis a ranger le resuitat dans la memoire 

20 m. Par convention, une telle operation sera notee de la fa$on suivante : 

m <— m + i 
Dans Toperation ci-dessus, Tincrement i peut etre : 

- soit une constante 

— soit le contenu d'une memoire 
25 — soit une expression arithmetique 

Dans Toperation ci-dessus, 1'increment i peut etre : 

• soit une constante 

• soit le contenu d'une memoire 

t soit une expression arithmetique 
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1.4.3 Procede d'ajustement dans lie cas de ponderateurs 
harmoniques 

Dans le cas ou les coefficients {u a (fc)} sont les ponderateurs harmoniques 
presentes dans la section 1.1.4 et dans Tannexe, les valeurs /*, f% et /* des 
5 fonctions d'ajustement /£(), /*() et /*() relatives au noeud <p a peuvent etre 
caiculees de fagon efficace a I'aide du procede suivant base sur le codage de 
5 : 



0) Allouer les memoires de travail suivantes 



10 



f x fy f 2 



i x lV f z 



15 



1) Effectuer les initialisations suivantes (l'ordre n'a aucune impor- 
tance) : 



<rj = = a% = 0 



2) Determiner le nombre n a de satellites de A(a). Ce nombre peut 
etre obtenu a partir du champ nb-sat de A(a) 

20 3) Soit A(a) 1 'ensemble des satellites de Fatome A(a) auxquels on 

accede par le champ sat de cet atome. Pour chaque satellite k 6 A(a), 
repeter les operations du bloc suivant : 

{ 
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3.1) Determiner le nombre de satellites de -A(Jk). Ce nom- 
bre peut etre obtenu a partir du champ nb jsat de A(k) 

3.2) Determiner les coordonnees {vtiVkiVk) ^ u noeU( * $k 
de Tatome A(k). Ces coordonnees peuvent etre obtenues a 

5 partir du champ noeud de A(Jfe). 

3.3) A partir des informations contenues dans A(k), calculer 
ou determiner la vaieur du ponderateur /z(Jk) (voir section 
1.1.4 et annexe). 

3.4) Effectuer les operations suivantes (l'ordre n'a aucune 
10 ' importance). 

% *~ + ti k ) 

+ 

$ r = s y a = 0 

3.5) Soit A(fc) l'ensemble des satellites de l'atome A(k) 
auquel on accede par le champ sat de cet atome. Pour 
chaque satellite ft € A(fc), /? ^ a effectuer les operations du 
bloc suivant : 

{ 

20 3.5.2) Determiner les coordonnees actuelles 

(<p0)<p y p,<pp) du noyau (pp de Patome A(f3). Ces co- 
ordonnees peuvent etre obtenues a partir du champ 
noeud de A({3). 

3.5.3) Effectuer les 3 operations suivantes (Pordre n'a 
2 5 aucune importance) : 

} 

3.6) Effectuer les 3 operations suivantes (I'ordre n'a aucune 
30 importance) : 

£<-/£ + K*)x xv?) 
/£<-/J + /*(*)xK-»*x^) 
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4) Effectuer les operations suivantes : 

M a *— ji(a) X n a X n a + <j m 

Mr «- M a 



5) Soit C(a) I'ensemble des contraintes attachees a I'atome A(a) 
auquel on accede par le champ const de cet atome. Pour chaque con- 
trainte c £ C(a) repeter les operations du bloc suivant : 

{ 

5.1) Determiner le type de la contrainte c 
1 0 5.2) En fonction du type de la contrainte c, effectuer les deux 

operations suivantes : 

5.2.1) Extraire les informations specifiques de la con- 
trainte c 

5.2.2) En fonction des informations extraites en (5.2.1), 
15 calculer les valeurs des expressions suivantes (voir sec- 
tions 2.5 et 2.7 de l'annexe) : 



20 



-pxa 
1 1 



t py<* r*za 
i 1 1 i t 

Qr,Q y i a >Q' a 

5.3) Effectuer les operations suivantes (l'ordre n'a aucune 
importance) : 



25 



} 



Ml' 



fz+-fz+ err - on 



6) Effectuer les operations suivantes (l'ordre n'a aucune importance) 



30 



(/*(«) X n a X <r% - fZ)/M t 
(/,(«) x n a X <rv - fV)/M? 
fa - W«) X »a X < - MM? 



n 



A la fin de ce procede de calcul, les memoires de travail /*, f% et /* contiennent 
les valeurs des fonctions d'ajustement /*(), f%{) et f*() utilisees pour ajuster 
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les valeurs des coordonnees du noeud $ a : 




5 Le procede d'ajustement decrit ci-dessus ne doit etre applique qu'aux noeuds 
dont la position est non fixee, c'est-a-dire aux noeuds qui ne sont pas des 
noeuds de controle (voir section 1.1.3). 



1*4.4 Variant es 



Variante 1 



1 0 La mise en oeuvre du procede d'ajustement (base sur la methode DSI) presentee 
au paragraphe precedent dans le cadre de ponderateurs harmoniques peut 
faciiement etre adaptee dans le cas de ponderateurs quelconques. La structure 
globale du procede reste la meme, seule change Tincrementation des differentes 
memoires de travail allouees a l'etape (0); la fagon d'incrementer ces memoires 

1 5 de travail depend des ponderateurs choisis. 

Par exemple, si nous choisissons les ponderateurs {v a (k)} suivants, 




si a = k 
si a € A(fc) 
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alors on verifie facilement que les fonctions d'ajustement associees sont : 

/:(•••) -tL-[ E *s)/ia(«)i 

~ E M*0 ■ (Evg) /I A(fc)| 2 - /|A(fc)|} 

- Eur -on ] 

M-l E <A)M*)\ 

H€A(a) ' 

E *»(*)• (E tf)/|A(k)P-tf/|A(fc)|} 

fe€A(a) "l^*' 

EW-QH ] 

/*(«)'( E rf)/|A(a)| 



10 



v *€A(a) 7 

- E #»(*)•{( E v>*)/|A(fc)| 2 -^/|A(fc)|) 

- eat -an ] 



avec 



y* 



*r-J4+Ei? 
^r=M a +E7? fl 

15 M a = (*(*)+ £ {M*)/|A(*)| 2 } 

*6A(a) 

Le procede de calcul des valeurs /*, et /* prises par ces fonctions au 
noeud <p a est alors quasiment identique a celui decrit pour les ponderateurs 
harmoniques : 
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0) Allouer les memoires de travail suivantes : 



fx fy f* 



t x *y % z 



1) Effectuer les initialisations suivantes (Pordre n'a aucune impor- 
tance) : 



10 



2) Determiner le nombre n a de satellites de A(at), Ce nombre peut 
etre obtenu a partir du champ nb_sat de A(ot) 

3) Soit A(a) Pensemble des satellites de l'atome A(ot) auxquels on 
accede par le champ sat de cet atome. Pour chaque satellite k € A(a), 

1 5 repeter les operations du bloc suivant : 

{ 
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3.1) Determiner le nombre nk de satellites de A(k). Ce nom- 
bre peut etre obtenu a partir du champ nb_sat de A(k). 
Effectuer ensuite Toperation suivante : 

n l *- n k X n k 



10 



15 



3.2) Determiner les coordonnees (p* du noeud (p k 
de l'atome A(k). Ces coordonnees peuvent etre obtenues a 
partir du champ noeud de A(k). 

3.3) A partir des informations contenues dans A(h) t calculer 
ou determiner la valeur du ponderateur /i(ib) (voir section 
1.1.4 et annexe). 

3.4) Effectuer les operations suivantes (1'ordre n'a aucune 
importance). 



$* = ,y= s *=o 



20 



25 



3.5) Soit A(Jfc) 1'ensemble des satellites de 1'atome A(k) 
auquel on accede par le champ sat de cet atome. Pour 
chaque satellite ft 6 A(Jb), (3 ^ a effectuer les operations du 
bloc suivant : 

{ 

3.5.2) Determiner les coordonnees actuelles 
{VfiiVfytVff) ^ u n °y au 9p & e I'atome A{fi). Ces co- 
ordonnees peuvent etre obtenues a partir du champ 
noeud de A(fi). 

3.5.3) Effectuer les 3 operations suivantes (l'ordre n'a 
aucune importance) : 



30 



} 



'%+v% 



3.6) Effectuer les 3 operations suivantes (1'ordre n'a aucune 
importance) : 



35 



rs <- n + m x w»i - *u»k) 
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} 



4) Effectuer les operations suivantes 



M a 



/*(<*) + * u 
M a 



5) Soit C(a) l'ensemble des contraintes attachees a I'atome A(a) 
auquel on accede par le champ const de cet atome. Pour chaque con- 
trainte c S C(a) repeter les operations du bloc suivant : 

{ 

10 5.1) Determiner le type de la contrainte c 

5.2) En fonction du type de la contrainte c, effectuer les deux 
operations suivantes : 

5.2.1) Extraire les informations specifiques de ia con- 
trainte c 

15 5.2.2) En fonction des informations extraites en (5.2.1), 

calculer les valeurs des expressions suivantes (voir sec- 
tions 2.5 et 2.7 de Tannexe) : 



20 



r>ra 



) 1 1 J i 

5.3) Effectuer les operations suivantes (l'ordre n'a aucune 
importance) : 



25 



} 



Mr 



M? + yf a 



Ml' 



M? + yf a 



fl + (r? a - QD 
£<-/S + (rr-Q?°) 



f y 

J a 



6) Effectuer les operations suivantes (l'ordre n'a aucune importance) 



30 



/*«- (/»(«) x<r5K-/*)/Af; 

/v^(Ma)x«rV/n 0 -/v)/MJ« 
K <- (,«(«) X **/n a - MM? 
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Les parametres numeriques stockes dans le champ info des atomes peuvent 
egalement etre interpoles par la methode DSI (voir section 2.1.2 et 2.2 de 
r annexe). Dans ce cas, la structure du precede d'interpoiation est quasi- 
ment identique a celle du precede d'ajustement de la surface presente ci- 
5 dessus; le seul changement consiste a ne conserver qu'une des trois fonctions 
d'ajustement normalement dediees aux trois coordonnees des noeuds et a cal- 
culer celle-ci en fonction du parametre interpoie. 

La presente invention trouve application tout 
particuli&rement dans les domaines suivants : 

10 - En g6ologie et en g6ophysique, elle permet de 

cr6er un module de surface, par exemple de const ruire une 
maquette solide en trois dimensions, sur laquelle on peut 
effectuer une simulation g6ophysique ou une simulation de 
1 'exploitation de richesses souterraines ; elle permet 

15 6galement d'obtenir un ensemble de donn£es propres <k 
initialiser un simulateur num6rique ou analogique de 
ph§nom6nes g6ophysiques ou d* exploitation de richesses; en 
outre, elle permet, & I'aide d'une station de travail 
graphique, d'afficher sur un 6cran ou d 1 imprinter sur un 

20 support des images permettant de contrdler les simulations 
ci-dessus; 

-en biologie ou en m6decine, 1* invention permet de 
construire des modules en trois dimensions d'organes et/ou 
de proth£ses, ou encore de simuler les effets d' operations 

25 en chirurgie plastique; de m§me, des images en deux 
dimensions affich6es ou imprim6es permettent ici encore de 
contrdler les operations ci-dessus. 

Mais il est clair que 1' invention s 1 applique plus 
g6n6ralement & la mod61isation de surfaces , de type 

30 quelconque, et notamment de toute surface naturelle ou 
industrielle* Elle peut trouver un int6r£t tout 
particulier dans le cadre de la conception assistee par 
ordinateur (CAO) . 

Enfin, la pr6sente invention n*est nullement 

35 iimit6e k la forme de realisation d6crite ci-dessus et 
r pr6sent6 sur les dessins, aais l'homme de l'art saura y 
apporter toute variante ou modification conforme k son 
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ANNEXE 



2. NOUVELLE VERSION DE LA METHODE DSI ADAPTEE A LA 

MODELISATION INTERACTIVE DES SURFACES COMPLEXES . 

Dans les domaines de ia Geologie et de la Biologic il est frequent d'avoir a modeliser des 
surfaces complexes correspondent par exemple a des interfaces entre zones de nature 

5 ou de proprietes differentes. Les techniques classiques de moderation basees sur les 
interpolations de Bezier ou sur les fonctions Splines ( voir ret [9]) sont tres mal adaptees 
a Ia prise en compte de ces donnees heterogenes. Sous proposons ici une approche 
differente basee sur la methode DSI (Discrete Smooth Interpolation * voir ret [8]). Dans 
cette approche, les surfaces sont modelisees avec des facettes triangulares irregu/ieres 

1 0 dont les sommets doivent etre determines aRn de tenir compte du plus grand nombre 
de donnees heterogenes possible. 



15 2.1 Introduction : 

Presentation des problemes etudies 

2.1.1 Surface 5 definitions assoctees 

20 Corame le montre la Figure 1 , soit S une surface composee de facettes triangulaires ad- 
jacentes plongee dans Tespace euclidien 3D (0|x, y, *). Cette surface nest pas for cement 
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n'est pas forcement connexe ni fermee. Dans ce qui suit, nous utiliserons ies 
notations suivantes : 

G = G{S) = ensemble des cotes des triangles 
ft = ft (5) = ensemble des sommets des triangles 
jV = jft| = nombre de points de ft 

G est un graphe dont les noeuds sont les points de ft et nous identifierons cet 
ensemble avec TenSemble des indices {l,-**,iV} utilises pour les numeroter. 
Le voisinage N(k) d'un noeud k € ft est defini par : 

N(k) = sous ensemble de ft constitue par les noeuds de 
10 ft qui peuvent etre atteints a partir du noeud k 

en s(k) pas au plus ou s(k) est un nombre positif 
donne. 

Dans ce qui suit, nous supposerons que ces voisinages N(k) possedent la pro- 
priety de symetrie suivante : 

15 a€iV(j3) <=> 0£N{a) 

Les voisinages ainsi definis engendrent une topologie sur G qui peut etre con- 
sideree comme une approximation de la topologie de S elle meme. 

2.1*2 Probleme 1 : Estimation <Tune fonction <p definie 
sur Q 

20 Supposons que la position de chaque noeud k 6 ft est connue et soit ip(k) 
une fonction definie pour chaque noeud € ft et connue seulement sur un 
sous-ensemble de ft : 

L = ensemble des noeuds / € ft ou = ip t est connu 

I = ensemble des noeuds i € ft ou <p(i) = <^ est inconnu 
25 = ft — L 

Generalement, L est different de ft et il y a une infinite de fonctions <p{k) 
definies sur ft et interpoiant les valeurs {<p{l) = <f>\ : I € L}. Notre objectif est 
de selectionner parmi toutes ces fonctions ceUe qui minimise un critere donne 
R*{ip) mesurant : 

30 • la rugosite globale de la fonction 

• Tecart entre la fonction et un certain nombre de contraintes et 
donnees qu'elle doit satisfaire 
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2.1*3 Probleme 2 : Ajustement de la surface S 

Dans ce cas, seuls les points cTun sous ensemble L de noeuds k € Q ont une 
position connue <p(k) = (p k : 

L = ensemble de noeuds / € Q ou = (p { est connu 

/ = ensemble des noeuds i € ft ou (p{i) = <£• est inconnu 
= Q-L 

Comme dans le cas du premier probleme, L est generalement different de Q, 
et il y a une infinite de surfaces S interpolant les points {£(/) = <p l : / € L}. 
Notre objectif est de selectionner parmi toutes ces surfaces, celle qui minimise 
1 o un critere #*(£) donne mesurant : 

• la rugosite globale de la surface S 

• l'ecart entre la surface S et un certain nombre de contraintes et donnees 
qu'elle doit satisfaire 

Compte tenu de Tindependance des composantes (p r (fc)> tp y {k), <p*(k)) de chaque 
15 vecteur ${k) y il est facile d'utiliser la solution du premier probleme pour 
resoudre le second en posant simplement : 

JF(£) = + + 

Dans le cas tres particulier ou i'ensemble Q correspond aux noeuds d'une 
grille rectangulaire posee sur la surface, on peut songer a utiliser une methode 
20 ^'interpolation du type spline, Briggs ou Bezier (Voir ref. [9], [1] et [3]); 
malheureusement ces methodes sont inappliquabies dans le cas des grilles 
irregulieres et ne permettent pas de tenir compte des contraintes presentees 
aux paragraphes 2.2.3, 2.5 et 2.7.2. 



2.2 Interpolation 

d'une fonction <p(k) definie sur ft 

Dans cette section, nous introduisons brievement la methode DSI (Discrete 
Smooth Interpolation) presentee dans l'articie [8J. 
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2.2.1 Notations 

Dans ce qui suit, nous noterons une matrice colonne de taille N telle que : 

Notre probleme ne depend pas de la methode de numerotation des noeuds de 
la grille. Afra de simplifier les notations, nous supposerons que les elements 
de la matrice <p ont ete eventuellement permutes de telle facon que <p puisse 
etre partagee en deux sous-matrices <pj et Vi telles que : 



10 



avec 



6/Va 



avec //j € i V |8 



!5 De plus, nous noterons || • ||£> une semi-norme associee a une matrice carree 
semi-definie positive [D] de taille (N x N) de telle facon que, pour toute 
matrice colonne X de taille N on ait : 

||* IB, = X*.[D]-X 

2.2.2 Definition d'un critere de rugosite globale R(<p) 

20 Soit R(tfi\k) le critere de rugosite locale defini par la formule suivante ou les 
{v°(fc)} sont des coefficients de ponderation positifs, negatifs ou nuls donnes : 

R(<p\k) =1 E vU 2 

R{<p\k) peut etre utilise pour construire un critere de rugosite giobale R{tp) 
defini de la fagon suivante ou fi(k) est une fonction de ponderation donnee 
2 5 definie sur ft : 

!>(*)•%>!*) 

k 

R{<p) est completement defini par la donnee des coefficients {^(fc)} et (i(k) et 
il est facile de verifier (voir ref. [8]) qu'il existe toujours une matrice (N x N) 
symetrique semi-definie positive [W] telle que : 



30 
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2.2.3 Prise en compte de contraintes lineaires au sens 
des moindres carres 

Considerons une matrice carree [A { ] de taille (N x N) et une matrice colonne 
B i de taille N definissant la contrainte lineaire suivante : 

5 [A { ]-<p ~ J9 { 

Etant donnee une matrice semi-definie positive [£>,*], le symboie figurant 
dans la formule ci-dessus a la signification suivante : 

II [A] — B{ Witi auss * P et i fc <l ue possible 

S'il existe piusieurs conditions de ce type, nous proposons de mesurer le degre 
1 o de violation de ces contraintes par le critere p(<p) defini comme suit : 

Mm>) = E II K-l 



2.2.4 Solution du probleme 

Parmi toutes les fonctions tp(k) definies sur ft et interpolant les donnees 
{<p(l) — ifii : / € L}, nous en choisirons une qui minimise le critere suivant : 

15 RW) = R(<P) + p(v) 

En developpant on obtient : 

R*(<p) = <p''[W*]-<p - 2'<p*.[Q] + C 
[W*\ = Z [ • [A] ' ML + [W] 

t 

20 C = Z II B i \\l 

i 

La partition de la matrice <p induit une partition similaire des matrices [W*} 
et Q utilisees pour definir : 



V = 



; Q = 



Qi 
Ql 
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La condition dR*(<p)/d(fi l = [0] implique l'equation suivante appelee "equation 
DSF et caracterisant toutes les fonctions {<p(k) : k € fl} solution de notte 
probleme : 

5 (DSI): < 

k avec : 0, = -[Wf L ] • y> L + Q, 

2.3 Unicite de la solution de Tequation DSI 
Definition 

L'ensemble L des noeuds ou tp est connu sera dit consistant relativement a 
1 0 R*(<p) si chaque composante connexe du graphe G contient au moins un noeud 
appartenant a L. 



15 



theoreme 

Si le critere de rugosite globale R(<p) est tel que 



1) L est consistant relativement a R(<p) 

2) (i(k) >0 VJb 

3) v a {k) >0 V {a ? k,<x £ N(k)} 

4) v k (k) < - 2 v tt (k) £ 0 



alors, Pequation DSI basee sur R((p) admet une solution unique. 



D emonst ration 

20 Le theoreme ci-dessus ne presente qu'un interet theorique et n'a aucune in- 
cidence sur le procede decrit dans le present brevet, c'est pourquoi nous ne 
donnerons pas ici sa demonstration. 

2 A Forme locale de l'equation DSI 

Dans ce qui suit, plutot que de resoudre directement Pequation DSI, nous 
25 proposons une approche iterative originale permettant d'eviter le caicui et le 
stockage de [W*]. Cette approche iterative est precisement celle qui est utilisee 
dans le present brevet. 
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Afin de simplifier les notations, nous supposerons dans ce qui suit que les 
matrices semi-definies positives [D { ] intervenant dans Tequation DSI sont 
diagonales. De plus, lorsqu'il y a des contraintes lineaires du type [A^ -ip ~ 
B i} nous utiliserons les notations suivantes : 

5 A* h = element de [A { ] situe sur la j ^^^ ne ^ 

m Qf = a th element de la matrice colonne ([A,-]' [£,-]£,•) 

2.4.1 Calcul de dW(v)/dv a 

Compte tenu de la definition de R(tp \k) y on peut ecrire Timplication suivante : 

m*) = E E W7 

/J€^(*)-r6iV(fc) 

10 = E (AW 2 -*! + E E 



dR(ip\k) 



' 2 • (v°(k)) 2 • Va + 2 • v«(fc) • 2 ' V/j si q € N(k) 



= < /»€«(*) 

i 0 sinon 



On deduit de cette relation que : 



1 dR(y) l v 3R(<p\k) _ 1 v ggMfe) 
15 2 2 A ^ a 2 fc€;v(o) «V a 

= V.- E *(*)(*•(*))' + E {?(*)""(*)' £ 
*€W(«) fc€N(<») w(*> 

Compte tenu de la structure diagonale de [DJ, il est facile de verifier que le 
a eme element de 5 1| [A,] ||2>. /a^ est tel que : 

20 = <?« • E A* fc (4 a ) 2 + E{ A fa A* • E 4%} - 0? 
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Compte tenu de la definition de on en deduit que 

1 



fy 0 



= ( p a 



{ E #)KW) 2 +EE^(^ S ) J } 
+ E to*) «■(*)• E A*)-^} 



*€AT(«) 



+ EEk^-EA^} 



- E<3f 

t 

La solution tp, est telle que dB?{<p)ldtf> a = 0» done la a eme composante yj a 
de v doit satisfaire l'equation OSI(a) suivante : 



10 



DSI(o) 



{ E W*) *■(*)• E At).^} 
+ Etff-Gf}} 



avec : < 



Ma= E A*) (*"(*))* 



15 Nous proposons d'appeler cette equation /orme /oca/e de l'equation DSI au 
noeud a. 



2.4.2 Proposition pour un algorithme iteratif 

La forme locale de l'equation DSI presentee dans la section precedente suggere 
un algorithme evident pour estimer la solution <p. Par exemple, a l'iteration 
20 la a eme composante <p£ +l ) de la solution <^ B+1 ' doit satisfaire l'equation 

DSI(a) ce qui permet de proposer le procede iteratif suivant : 



s it 7 l'ensemble des indices de noeuds ou <p a est inconnu 
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8 it <p une solution approchee initiale 

tant_que( il est necessaire de faite de nouvelles iterations ) 

{ 

pOUT_tOUt( Of e / ) 

{ 

Ml = {•••} 

5 = •{'•'} 

}• 

} 

Comme on peut le voir, cet algorithme tres simple n 'utilise pas explicitement la 
matrice [W*] intervenant dans Tequation DSI mais il requiert de recalculer un 
grand nombre de fois des produits du type {v a (k)*v^(k)} qui sont precisernent 
ceux utilises pour construire [W]. 

1 0 En fait, si la solution approchee initiale est proche de la solution exacte, 
alors il n'est pas necessaire d'effectuer beaucoup d'iterations et ie surcout de 
calcul devient negligeable. Par exemple, cette situation se presente dans un 
contexte interactif ou la solution initiale <p est prise egale a la solution avant 
que quelques modifications locales soient effectuees par Tutilisateur. Dans ce 

1 5 cas, malgre le leger surcout engendre, la forme locale de Pequation DSI est 
preferable car elle est bien plus facile a mettre en oeuvre que la forme globale. 



Afin de montrer que Talgorithme propose converge, on peut remarquer 
que peut etre exprime sous forme d'une fonction de ip a ayant la forme 

suivante : 

0 IT{tp) = A-tpl + B-<p a + C 

Les coefficients A, B et C sont independants de ip a et on a : 

b = 2. { s w*) *■(*)■ e ^w.^+Ew-on] 

Le minimum du critere global IFfa) est atteint pour la valeur tp a = -B/(2A) 
25 qui correspond precisernent a la valeur fournie pax l'equation DSI(a). Ceci 
nous permet de conclure que notre algorithme converge car, a chaque iteration 
du processus iteratif, la valeur de la fonction positive ou nulle R*(ip) decroit : 
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2.5 Prise en compte de donnees imprecises 

Dans ce qui suit, nous introduisons deux types de contraintes lineaires qui sont 
d'une extreme importance en moderation geometrique ainsi qu'on le verra 
au paragraphe 2.7.6. L'objectif de cette section est de prendre en compte des 
5 donnees imprecises concernant les valeurs inconmies {(p\ : A € /}. De telles * 
donnees sont supposees associees a des coefficients de certitude & 2 6 R + qui 
sont en fait des coefficients de ponderation positifs. 



2.5*1 Cas d'une donnee ponctuelle imprecise 

Soit A € / un indice de noeud pour lequel la valeur inconnue <p x est supposee 
1 0 proche d'une donnee incertaine <£ A avec un coefficient de certitude donne egal 
a xs x : 

ou ~ correspond a la condition suivante : <p x : 

xsj\ • \tp* — $ x \ 2 minimum 

1 5 Une telle condition est aisement prise en compte par la methode DSI sous 
forme d'une i eme contrainte avec : 

+1 si fc =s h = A _ f <p x si k = A 

0 sinon * "I 0 sinon 




20 , v smon 



n kh _ f *\ si 
U{ " \ 0 sir 



. „ M , ~ , - k =s h = A 
et 



Les coefficients Tf , et 7? correspondants utilises dans i'equation locale 
DSI(a) sont alors definis par : 



17 = 0 Qf = { n 
1 1 [ 0 smon 



w\ • (p x si a = A 



w\ si a = A 
et 7T = i o* sinon 



et 7 ? = { 
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2.5.2 Cas d'une donnee differentidle imprecise 

Soit A € / un index de noeud pour lequel ia valeur inconnue <p x est supposee 
liee a une autre valeur correspondant a un autre index fi ^ A, Supposons 
que ce lien ait la forme suivante ou A Xfi est une valeur donnee : 

5 (Vp - V\) - ^a m 

Cette relation est supposee imprecise et est interpretee de la fa?on suivante 
ou est un coefficient de certitude : 

wJn ' Krt - <P\) - a a m I 2 minimum 

Une telle condition peut etre facilement prise en compte par la methode DSI 
10 en tant que i eme contrainte avec : 



15 



' -1 


si k = 


A = A 




-1 


si k = 


h = fi 




+1 


si k = 


A et h 


= t 1 


+1 


si k = 


fi et h 


= A 


0 


sinon 







[ + A Am si * = A 
- si k = 

0 sinon 







si k = 


et Dj h = < 


1—2 


si fc = 




0 


sinon 



Les coefficients T", (J* et yf correspondants utilises dans l'equation locale 
20 DSI(a) sont alors definis par : 



( -*>ln"t>l sia = A 
r? = i -^t * ¥>a si a = p. 

0 sinon 



-^L'Aam sia = A 



0 



sinon 



f _2 
C7 



25 



et 7 f=< 



A „ sia = A 
»L sia = p 
0 sinon 
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2.5.3 Choix des coefficients de certitude 

Considerons le terme M* de I'equation DSI(a) : 

M; = M ff + 7f + 7 f + ••• 
avec:M 0 = £ p(k) (v a (k)) 2 

Dans les deux exemples presentes ci-dessus, les termes ff sont soit egaux a 
zero, soit correspondent au coefficient de certitude : 

7f = *1 

Cette remarque nous suggere de choisir pour us\ un pourcentage donne p- de 

Le terme Af* de I'equation DSI(or) est alors egal a : 

K = M a '(l+ Pl {a) + p 2 {a) + >") 
ou Pi(a) est soit egal a un pourcentage donne p t - ou egal a zero. 



is 2*6 Choix des coefficients de ponderation 

Le choix des coefficients de ponderation {v a (fc)} et {fi{k)} est compietement 
libre excepte que les coefficients {fi{k)} doivent etre positifs ou egaux a zero. 
Dans ce qui suit nous donnons un exempie de choix de ces deux families de 
coefficients. 



20 2.6.1 Choix de coefficients {v Q (k)} : 
ponderation harmonique 

Soit A(fc) T'orbite" de k definie de la fa$on suivante : 

A(fc) = N(k) - {k} 

Soit |A(A:)| le nombre d'elements de A(fc). Nous dirons que nous utilisons une 
25 ponderation harmonique si les coefficients {v a (k)} sont choisis de la fa^on 
suivante : 

ti*r«- J H A (*)I si a = fc 
v W " \ 1 si a e A(Jb) 
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Les fonctions harmoniques ont ia propriete caracteristique d'etre egales en tout 
point a leur propre moyenne sur un cercle centre en ce point. Comme on peut 
le voir, R{<p\k) s= 0 si tp k est egal a ia moyenne des valeurs de (p a entourant 
le noeud k, c'est pourquoi, par analogie avec les fonctions harmoniques, nous 
5 proposons d'appeller harmoniques les coefficients {^"(A:)} presentes dans cette 
section. 



Si nous developpons l'equation DSI(a) avec ces coefficients, nous obtenons 
la formule suivante qui peut etre facilement transcrite dans un langage de 
programmation : 



10 



DSI(a) 



/«W-|A(a)|.(E <P>) 
- EiTf-Qf} } 



avec 



*€A(a) « 



15 2.6.2 Choix des coefficients {/*(*)} 

Les coefficients {p(fc)} ont ete introduits afin de moduler localement le lis- 
sage de Interpolation. En Pabsence de raison particuliere, nous suggerons 
d'utiliser une 
ponderation unifonne : 

20 Mfc) = lVfc€fi 

Parmi toutes les ponderations non uniformes, la ponderation suivante semble 
particuiierement interessante s'ii est necessaire devoir une fonction lisse au 
voisinage des donnees : 



25 




1 si k € / 
m > 1 si k€ L 
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2.7 Application a la 

Modelisation Geometrique 



Dans cette section, nous nous proposons de montrer comment utiiiser la 
methode DSI afin d'estimer la surface trianguiee S elle meme pour laqueile 
5 on suppose disposer des informations partielles suivantes : 

# la position exacte de quelques sommets de triangles, 

• la position approchee de quelques sommets de triangles, 
« des contraintes vectorielles pr^cisant la forme de «S, 



1 o Pour cela, on est amene a definir ip a comme le vecteur courant joignant 
Torigine de R 3 au noeud courant de SI et nous noterons <p x , <p y et ip z ses 
composantes sur la base orthonormee {ox,oy,o£} de R 3 : 



2*7.1 Definition (Tun critere de rugosite globale R(<p) 

1 5 Suivant notre habitude, Pensemble Q, des noeuds de la grille est partage en 
deux sous-ensembles I et L 



Les noeuds J 6 L dont la position $ 6 R 3 est connue seront appeles noeud 
20 de controle et Ton se propose de determiner la position des points restants 
i € / de telle fagon que le critere local i?(<p|fc) defini ci-dessous et base sur un 
ensemble donne de ponder ateurs {v a (k)} soit aussi petit que possible : 



Ce critere constitue la forme vectorielle du critere DSI local. Le critere DSI 
25 global associe R($) est alors defini de la fagon suivante ou les {(i{k)} sont des 
coefficients de ponderation non negatifs donnes : 



• etc . . . 



(p { est inconnu V i € / 
9i est connu V I € L 



k 
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Si l'on definit R(<p x ), ^(v 9 ) e t ^(v 2 ) a i ns i qu'il aete fait dans les paragraphes 
precedents : 

I W) = !>(*)• avec :j R(^|fc) = | £ «•(*). p 

fc «6/V(Jfe) 

*(v> 2 ) = E^W'W) avec : %>*|fc) = | X) t/*(*)-? 2 ! 2 
alors, en utilisant le theoreme de Pythagore, on verifie faciiement que Ton a : 



2.7*2 Prise en compte de contraintes lineaires au sens 
des moindres carres 

1 0 Considerons les matrices suivantes dimensionnees de telle f agon que N represente 
le nombre total de noeuds de Pensemble fl : 



Vi/ = x,y y z : < 



[D u iv ] . est une matrice diagonale non negative de taille N X N 
[A v iu ] est une matrice carree de taille N X N 



Bi est une matrice colonne de taille N 



1 5 Si nous voulons que les composantes <p z y <p y et <p 2 satisf assent les contraintes 
suivantes 

[Att'9' * Bl 
20 alors on est amene a introduire le critere suivant 

p{$) = p{ V *) + ,,(v> v ) + P(<fi') 

f = E IIKl-^-uEIIlt 

avec:< - £ II W ^ " ^ »H 

rtv 2 ) = E ll^j-v 2 -* 2 ,!!^ 
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2.7.3 Solution du probleme 

Parmi toutes les surfaces S interpolant les donnees {<£(/) = <p t : / € L), nous 
nous proposons de selectionner celle qui minimise le critere suivant : 



5 II est facile de verifier que Ton peut toujours decomposer cette formule de la 
fagon suivante : 



Compte tenu de Independence des composantes <p x , (p y et <^*, on en conclut 



La minimisation du critere vectoriel DSI est done equivalente a la minimisation 
des trois criteres DSI correspondant aux trois composantes du vecteur courant 

2.7.4 Remarque 

20 Si la surface S est composee de plusieurs morceaux disjoints, alors le theoreme 
present e au paxagraphe 2.3 nous garantit que notre probleme admet une so- 
lution unique des iors qu'il y a au moins un sommet de facette ayant une 
position connue sur chaque morceau. 

2.7.5 Utilisation interactive 

25 En modelisation geometrique, Tobjectif est de definir interactivement la po- 
sition dans Tespace du graphe G associe a Tensembie des noeuds SI Pour ce 
faire, on suppose connue une forme initiate et Tutilisateur deplace quelques 
noeuds et modifie interactivement un certain nombre de contraintes (noeuds 
de controles et donnees imprecises); ces modifications sont alors propagees 

30 aux noeuds i € I en utilisant la forme vectorielle de la methode DSI. 



( RT{f) = R{<p*) + p(v*) 
avec: I = R(tp y ) + p(<p y ) 



que : 



f minimum 
/T(9) minimum I R?(<P y ) minimum 
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Nous avons montre dans les paragraphes precedents que ce probieme se 
decomposait en trois sous-problemes correspondant aux trois coordonnees 
(x,y,z) et la meilleure fagon de resoudre ces trois problemes est d'utiliser 
Talgorithnie iteratif associe a la forme locale des equations DSI correspon- 



D'une etape a Tautre du processus de modelisation on peut modifier non 
seulement la position des coordonnees des noeuds de controle mais Tensemble 
L lui meme peut etre change si necessaire; par exemple, si nous voulons appli- 
quer DSI seulement a un sous ensemble T de noeuds, alors nous devons choisir 
10 L de telle fa^on que L D f ou f designe Tensemble complementaire de T dans 



De plus, compte tenu de la decomposition de notre probieme en trois sous- 
problemes resolus independemment et correspondant aux trois coordonnees, 
on peut, si necessaire, definir des sous ensembles differents L x > L y et L 2 pour 
1 5 chacune de ces coordonnees; par exemple, si les coordonnees (2, y) ne doivent 
pas etre modifiees, alors il suffit de prendre L x et L y egaux a Tensembie Q 
tout entier. 

• 2*7.6 Prise en compte de donnees geometriques impr^cises 

Initialement, la methode DSI a ete con§ue pour modeliser les surfaces com- 
20 plexes rencontrees dans le domaine des sciences naturelles comme par exemple 
la geologie et la biologie. Dans ce cas, on doit generalement prendre en compte 
de nombreuses donnees impresices comme pax exemple : 



25 Noeuds de controle flous 

Par definition, nous appellerons "noeud de controle flou" tout noeud A 6 
I dont la position tp\ est connue comme etant approximativement egale a 
un vecteur donne (p^ et ceci avec un certain degre de certitude vo\. En se 
referant au paragraphe 2.5.1, on voit que Ton peut prendre en compte une 
30 telle information imprecise a Taide de la methode DSI en introduisant les 
contraintes suivantes dans les criteres R?^) et ^ t { i P 2 ) : 



5 dantes. 



{ 



• des noeuds de controle flous 

• des donnees vectorielles floues 
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Contraintes vectorielles floues 

Dans de nombreuses applications, il est necessaire de controler le vecteur 
(<? M — <? A ) joignant deux noeuds <£ A et tp^ d'une surface. Par definition, une telle 
contrainte sera appellee "contrainte vectorielle floue w si elle doit etre satisfaite 
5 avec un certain degre de certitude gt Am . Quelques exemples de contraintes de 
ce type sont presentes sur la Figure 2. 

En pratique, au moins un des deux points £ A ou (p^ a une position inconnue, 
par exemple celui correspondant a A € J, et parmi les cas importants, on peut 
citer les deux suivants : 

10 • Le premier cas correspond aux Figures 2. a et 2.b ou Ton desire controler 
la forme d'un objet de telle fagon que (<p M - <p x ) soit egal a un vecteur 
donne A A|i . 

• Le second cas correspond a la Figure 2.c ou I'on desire controler la forme 
d'un objet de telle fagon que chaque (jPn ~~ $x) so ^ orthogonal a un 
1 5 vecteur donne V* Cette situation se presente lorsque Ton veut controler 

le vecteur normal a une surface; dans ce cas, si V est le vecteur normal 
au noeud A de la surface, il est judicieux de controler son orthogonaiite 
avec chacun de ces vecteurs (<£ M - £ A ) pour tout p 6 AT(A). 

Si Ton se reporte au paragraphe 2.5.2, le premier cas est facilement pris en 
20 compte par la methode DSI en introduisant les contraintes suivantes dans les 
criteres -R*(v y ) et avec un coefficient de ponderation egal a 

Pour le second cas, le probleme n'est pas tres different car on a : 

W - vl) * -{v* • - v5) + V' • (<p'„ - <pi)}/v* 
M - vl) * -{V s • (9* - + v» • (<p* - vl)}/v* 

Dans ce cas, nous suggerons de prendre les valeurs des termes (va»Va>vI) et 
(Vm»Vm>Vm) i^tervenant dans les membres de droite de ces equations egales 
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aux valeurs correspondantes obtenues a l'iteration precedente; pour se ramener 
au cas precedent, il suffit alors de poser : 

= -{v»-M-vS) + v».(y;-9j)}/v 
5 a* Am = -<v.M-v«) + v».w-v5)}/v 

Remarque 

Dans la prise en compte des donnees imprecises presentee ci-dessus, nous 
n'avons utilise qu'un seul coefficient de certitude ts^ M commun aux trois com- 
posantes de ces contraintes. En pratique, ceci n'est pas une obligation et, 
1 o compte tenu de I'independance des trois composantes, il est toujours possible 
d'utiliser trois coefficients de certitude differents G7*£, cct^ et cj^ correspon- 
dant a ces trois composantes. 

2.8 Generalisations 

Avant de conciure, mentionnons quelques generalisations evidentes : 

15 • La methode DSI peut etre utilisee avec n'importe quels types de facettes 
polygonales. En pratique, nous conseilions d'utiliser des triangles car ce 
sont les polygones les plus simples et les plus faciles a manipuler dans 
un ordinateur. 

• La methode DSI peut etre utilisee pour estimer des courbes polygonales 
20 dans un espace 3D; dans ce cas, les triangles sont remplaces par des 

segments. 

2.9 CONCLUSION 

La forme locale de Tequation DSI constitue un outil simple et puissant pour 
modeliser les surfaces complexes que Ton peut rencontrer par exemple en 
25 geologie ou en biologie. En particulier : 

• Le maillage utilise pour modeliser une surface peut etre compietement 
libre et il est possible d'utiliser des algorithmes automatiques pour les 
construire; par exemple, si le maillage est compose de triangles, on peut 
utiliser des algorithmes derives de la methode de Delaunay (voir ref. 

30 [2] et [5]). De plus, la taille du maillage peut facilement etre ajustee 

localement en fonction de la complexity de la surface. 
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La methode DSI permet de prendre en compte un tres grand nombre 
de donnees precises ou imprecises. Cette possibility est particulierement 
interessante dans le domaine des sciences natureUes ou I'objectif n'est 
pas de generer des surfaces esthetiques mais d'ajuster ces surfaces aux 
5 donnees precises et imprecises de fagon coherente et efficace. 

• I/aigorithme derive de Tequation DSI locale est suffisamment rapide 
pour permet tre une utilisation interactive. 

Nous n'avons rien dit sur la representation des facettes elles-memes, Selon ies 
applications, on peut prendre des facettes planes mais, si cela est necessaire, il 
1 0 est toujours possible de representer ces facettes par des morceaux de surfaces 
non planes interpolant les noeuds du maiilage (voir ref. [4]). 
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REVENDICATIONS 

1. Proc6d6 d'obtention d'une mod61isation d'une 
surface (S) repr6sentative notamment de 1' inter face entre 
deux zones de natures ou de propri6t6s diff6rentes dans 
un corps & trois dimensions tel qu'une formation 
g6ologique ou un corps vivant, caract6ris6 en ce qu'il 
comprend les Stapes consistant k: 

- obtenir au moyen d'appareils de mesure un 
ensemble de donn6es g6om6triques relatives k la surface 
et associ6es respect ivement k des points donn&s de ladite 
surface? 

- effectuer un maillage de la surface de telle 
sorte que l 1 ensemble desdits points donnas soit un sous- 
ensemble des noeuds (tj^ du maillage; 

- m6moriser k une adresse donn&e d'une m&moire, 

sp6cifique k chaque noeud { t ) du maillage, les 

informations suivantes: 

at U £ 

* les coordonn6es (Cj^, ^) du noeud consid6r6, 

* le nombre (nb-sat) de noeuds satellites du 
noeud consid6r6, 

* une information (sat) permettant d'acc6der aux 
20 adresses sp6cifiques desdits noeuds satellites et par 

suite aux informations relatives k ces derniers, 

* le cas 6ch6ant des donnges g6om6triques 
(const) associ6es audit noeud consid6r 6 , 

- pour chaque noeud du maillage, d6finir un indice 
25 de rugosit6 local (R(C^jk)) obtenu k partir d'une somme 

pond6r6e des coordonn6es actuelles du noeud et de ses 
satellites, 

- d6finir la somme (R* «^) ) d'un indice de rugosit6 
global (RKp)), constitu6 par 1' addition des indices de 

30 rugosity locaux associ6s k chaque noeud , et d'un indice 
global (^>(<f)) de violation desdites donn6es g6om6triques; 

- ajuster de fagon it6rative les coordonn6es de 
chacun des noeuds dont on ne connait pas les coordonn6es 
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precises, n utilisant a chaque fois 1' addition d'une 
combinaison ponderee des coordonne.es actuelles des 
satellites et des satellites des satellites dudit noeud 
et d'une combinaison des donn6es g6om6triques associfees 
audit noeud, de mani&re 4 amener ladite somme 4 un 
minimum; 



- cr6er une repr6sentation de la surface 4 partir 



2. Proc6d6 selon la revendication 1, caract6ris6 
en ce qu'on effectue un maillage triangulaire de la 
surface, 

3. Proc6d6 selon l'une des revendications 1 et 2, 
caract6ris6 en ce que l'6tape d'ajustement comprend trois 
sous-6tapes distinctes mettant en jeu respectivement les 
trois coordonn6es spatiales (<^, (fy,C^) des noeuds et des 
satellites, 

4. Proc6d6 selon l'une des revendications 1 a 3, 
caract6ris6 en ce que les donn6es g6om6triques 
comprennent des donn6es de coordonn6es (<fk ) des noeuds . 

5. Proc6d6 selon l'une des revendications 14 4, 
caract6ris6 en ce que les donq 6 es ^ g6om6triques 
comprennent des donn6es de vecteurs I Cj>^ - *f u ) entre deux 
noeuds donn6s. 

6. Proc6d6 selon l'une des revendications 14 5, 
caract6ris6 en ce que les donn6es g6om6triques 
comprennent des donn6es de vecteur normal (V) 4 la 
surface k mod61iser. 

7. Proc6d6 selon l'une des revendications 4 4 6, 
caract6ris6 en ce qu'au moins certaines desdites donn6es 
g6om6triques sont assorties d'un coefficient (^v^ty) 
repr6sentatif du degr6 de certitude avec lequel ces 
donn6es sont connues et en ce que l'indice de violation 
des donn6es intdgre lesdits coefficients de certitude, 

8. Proc6d6 selon l'une des revendications 
pr6c6dentes, caract6ris6 en ce qu'il comprend en outre 
les 6tapes consistant 4 mesurer des propri6t6s du corps 4 




chacun des noeuds. 
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trois dim nsions dans la region d'au moins certains 
noeuds du maillage t k m6moriser les mesures obtenues 
(info) aux adress s sp6cifiques aux noeuds consid6r£s. 

9. Proc6d£ selon l'une des revendications 
5 pr6c6dentes, caract£ris6 en ce que, dans chaque 

coefficient de rugosit6 local, le rapport entre le poids 
(v« (k) ; a=k) du noeud concern^ et le poids (v° (k) ; ae 
A.(k)) d'un satellite quelconque dudit noeud est 6gal au 
nombre (nk ) de satellites dudit noeud, affecte d'un signe 
1Q n6gatif. 

10. Proc6d6 selon la revendication 9 f caract6ris6 
en ce que ladite combinaison pond6r6e utilise les memes 
poids. 

11. Proc6d6 selon l'une quelconque des 
15 revendications pr6c£dentes, caract6ris6 en ce que la 

sorome des coefficients de rugosite locaux est une somme 
pond6r6e . 

12. Proc6d6 selon la revendication 11, rattach6e & 
la revendication 10, caract6ris£ en ce que ladite 

20 combinaison pond6r6e utilise 6galement les poids (^(k)) 
associ£s aux coefficients de rugosit6 locaux. 

13. Proc6d6 selon l'une des revendications 
pr6c6dentes,, caract6ris6 en ce que la surface est une 
surface der separation entre deux milieux de propria t£s ou 

25 de natures diff£rentes dans une formation g£ologique et 
en ce qu'on utilise la ra?r6sentation de la surface pour 
optimiser 1 1 exploration de richesses souterrames de 
ladite formation. 

14. Proc6d£ selon l'une des revendications L a 12, 
caract6ris6 en ce que la surface est representative des 

30 variations d'un propri£t6 d'un corps k deux dimensions au 
voisinage d'un plan traversant ledit corps. ladite 
propri6t6 6tant representee par la dimension restante, et 
en ce qu'on utilise la representation de la surface pour 
optimiser 1 'exploration ou 1' exploitation ' de richesses 

35 souterraines d ladite formation. 
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15. Procede s Ion l'une d s revendications 1 a 14, 
caract6rise en ce que la repres ntation de la surface est 
une repr6s ntation graphiqu sur un support plan. 

16. Proc6de selon l'une des revendications 1 4 14, 
5 caract6ris6 en ce que la representation de la surface est 

une maquette en trois dimensions. 

17. Dispositif pour I'obtention d'une mod61isation 
d'une surface fS) representative notaroroent de l'interf ace 
entre deux zones de natures ou de propri6t6s diffferentes 
dans un corps k trois dimensions tel qu'une formation 
geologique ou un corps d'etre vivant, caract£ris£ en ce 
qu'il comprend : 

- des moyens de mesure un ensemble de donn6es 
g6om6triques relatives k la surface et associ6es 
respect ivement k des points donn6s de ladite surface; 

- des moyens pour effectuer un maillage de la 
surface de telle sorte que 1' ensemble desdits points 
donn&s soit un sous-ensemble des noeuds du maillage; 

- une ra6moire contenant k des adresses donn£es r 
sp6cifiques k chaque noeud du maillage, les informations 
suivantes: 

* les coordonn6es (Cffc^ Cfjf#(f du noeud consid6r6, 

* le nombre (nb-sat) de noeuds satellites du 
noeud consider 6, 

* une information (sat) permettant d'acc6der aux 
adresses sp6cifiques desdits noeuds satellites et par 
suite aux informations relatives k ces derniers, 

* le cas 6ch6ant des donn6es geometriques 
(const) associ6es audit noeud consider, 

- des moyens de calcul pour ajuster de fagon 
30 iterative les coordonn6es de chacun des noeuds dont on ne 

connait pas les coordonn6es precises, en utilisant k 
chaque fois 1' addition d'une combinaison pond6r6e des 
coordonn6es actuelles des satellites et des satellites 
des satellites dudit noeud et d'une combinaison des 
35 donn6es g6om6triques associ6es audit noeud , de maniere k 
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amener 4 un minimum une aomm (RMCf)) d'un indice de 
rugosit6 global (R(Cf)), constitu6 par l 1 addition 
d'indices de rugosit6 locaux (R ( Cp | k ) ) respectivement 
associ6s k chacun des noeuds et chacun obtenu k partir 

5 d'une somme pond6r6e des coordonnSes actuelles dudit 
noeud et de ses satellites, et d'un indice global (^>(<f)) 
de violation desdites donn6es g6om6triques; et 

- des moyens pour cr6er une representation de la 
surface k partir des coordonn6es ajust6es de chacun des 

10 noeuds . 
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